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Abstract 

We re-visit the problem of New Physics (NP) contribution to the branching ratio of B s — >• 
decay in light of the recent observation of this decay by LHCb. We consider R-parity 
violating (RPV) supersymmetric models as a primary example - recently one has reported stringent 
constraints on the products of the RPV coupling constants that account for the B s -4 
transition at the tree level. We argue that despite the LHCb measurement of B(B S — > u + u~) is 
in a remarkable agreement with the Standard Model (SM) prediction, there is still a room for a 
significant New Physics contribution to B(B S — > fi + u~), as the sign of the B s — > transition 
amplitude may be opposite to that of the Standard Model, or alternatively the amplitude may 
have a large phase. We conduct our analysis mainly for the case of real RPV couplings. We find 
that taking into account the scenario with the sign flip of B s — > [J> + u~ amplitude (as compared to 
that of the SM) makes the bounds on the RPV coupling products significantly weaker. We also 
discuss briefly how our results are modified if the RPV couplings have large phases. In particular, 
we examine the dependence of the derived bounds on the phase of the NP amplitude. 
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The rare B s — > decay is believed to be one of the most powerful tools to test the 

physics that may occur beyond the Standard Model. Within the Standard Model this decay 
is loop-induced and in addition is helicity suppressed. Numerical evaluation gives [l-3] 

B(B a ->/x + aO = (3.25 ±0.17) x 10" 9 (1) 

In contrast, B s — > fi + fi~ decay rate may be dramatically enhanced within some of the 
Standard Model extensions and may exceed the SM prediction by several orders of mag- 
nitude. At the same time this decay is characterized by a pure final leptonic state, which 
causes the theoretical predictions for it to be very clean. It was therefore used intensively to 
constrain the SM extensions, and there was a hope to observe a distinct New Physics signal 
in this decay mode. 

Recently the LHCb collaboration has reported the first evidence for B s — > /U + yU~ decay 
at 3.5cr level |4j, 

B exp (B s -)■ fi+fi-) = (3.5±J;1) x 1(T 9 (2) 

which is in a remarkable agreement with the SM prediction. However, it would not be correct 
to declare that there is no New Physics contribution to B s — > fi + n~ at all. Some of the 
popular SM extensions do predict indeed a negligible NP contribution to B s — > n + n~ decay 
rate, due to strong correlations between the B s — B s mixing and B s — > fi + fi~~ amplitudes 
jsj]. Yet, for other SM extensions the problem of New Physics contribution to B s — > fi + n~ 



in light of t 



re recent observation of this decay by LHC 



} is the subject of discussion in the 



literature {^[gU^]. In particular, it has been argued in {6-8] that the LHCb result still leaves 
a room for a non-negligible NP contribution, due to the uncertainty in the experimental 
value of B(B S — > 

In this paper we examine a source of New Physics contribution to B s — > that would 

be actual even in the idealized limit of zero experimental and theoretical uncertainties in 
B(B S — > and perfect coincidence of the SM prediction with the experimental data. 

Namely, we consider a possibility for the B s — > fi + {i~ transition amplitude to have a sign 
opposite to that of the Standard Model or to have a large phase. The LHCb measurement 
of B s branching ratio constrains the decay rate, whereas the sign or phase of the 

transition amplitude remains arbitrary. Thus, it is possible that: 

• If the amplitude is real (or has a small enough phase so that it may be discarded), one 
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may fit the experimental data for B(B S — > ) in particular when 

A NP (B S ~ -2A SM (B S -> /i + /0, so that 

->• //+/0 = ^(5, -> + A SM (B S -> //+//-) ~ -A 5A/ (B S -> /iV) (3) 

• If instead the amplitude has a large phase, one may fit the experimental data for 
B(B S when 

|A SM (B a -> /iV) + l^ P (5 s -> /iV")|e^| ^ |A 5M (5 S -> ^V)| (4) 

(if neglecting the SM amplitude phase and using the approximation B exp (B s — > 
/ttV - ) ~ B SM (B S -> Note that Eq. (gj implies 

- 2A SM (5 S < i?e [^ P (5 S ->■ < 0, (5) 

\lm [A NP {B S Ai+A* - )] | < |^ 5M (5 S A*V)| (6) 

In particular, 

|Jm [A^ p (B a ->• | ~ |A 5M (5 S ->■ if 

i?e [A JVP ( J B S -> M V)] ~ -A SM (5 S -> (7) 

One may infer from Eq.'s ([3]) - (TTJ) that the NP contribution to the B s — > transition 
amplitude is the largest when the amplitude just flips the sign as compared to that of the 
Standard Model (rather than gets a large non-trivial phase). So, we will be concentrating 
here mainly on the case of a real amplitude, by assuming that the relevant NP parameters 
are real. We will however discuss at the end of the paper how our results are modified in 
presence of large phases of the NP parameters. 

Note that the possibility of the B s — > amplitude sign flip has already been men- 



tioned in {7] where one considered New Physics models with modified Z-boson couplings to 
down-type quarks. This possibility has been rejected there, as it is disfavored by the con- 
straints onZ-s- bb. To our best knowledge, there is no reason to disfavor the B s — > 
amplitude sign flip within other SM extensions (in fact it has also been implicitly considered 
in |2| within the general analysis of the NP contribution to B s — > ^ in the variety of the 
models with the amplitude phases varied freely form to tt). In our opinion, the detailed 
analysis of the possibility that B s — > amplitude may has a sign opposite to that of the 



Standard Model (and/or have a large phase) may be of great importance, especially in light 
of future improvement of the experimental accuracy of measurements of the B s — > fi + fi~ 
branching ratio. 

We consider here R-parity violating supersymmetric models with leptonic number viola- 



tion as a primary example. It has been recently argued |q| that the remarkable agreement 
between the LHCb measurement and the SM prediction for B s — > branching ratio 

implies rigorous constraints on the RPV coupling products that account for the B s — > fi + n~ 
transition at the tree level. We show that if the B s — > transition amplitude is allowed 

to have a sign opposite to that of the Standard Model, bounds on the RPV couplings may 
be by order of magnitude weaker. 

The most general Yukawa superpotential for an explicitly broken R-parity supersymmet- 
ric theory may be written as 

W M = -XijkULjEl + X' ijk LiQjDl + -\" jk U?D?Dl (8) 

Here Q and L denote SU (2) L doublet quark and lepton superfields, and U, D and E stand 
for the SU{2) L singlet up-quark, down-quark and charged lepton superfields. Also, i, j, k = 
1,2,3 are generation indices. We shall require baryon number symmetry by setting A" jfc to 
zero. Also, as mentioned above, we will assume the couplings A^-fc and A^- fc are real. 

Subsequently, the Lagrangian describing the RPV SUSY contribution to B s — > can 
be written as 

£r = - ^A- 23 Z^kP L S + A-32Z}i^ (9) 

where Pl,r are the helicity projection operators, and we use the notation — (1 — 7s)/2. 
Note that for the sake of transparency of our analysis, we neglect the transformation of the 
RPV couplings from the weak isospin basis to the (s) quark and sneutrino mass basis. (We 
invoke however the reader to be cautious when using the bounds on RPV couplings derived 
in this paper. Rigourously speaking, they may be used for the processes involving down 
type quark - down type quark - sneutrino and down type quark - up type squark - charged 
lepton transitions only.) 

Within R-parity violating supersymmetric models, to the lowest order in perturbation 
theory B s — > fi + fi~ transition occurs at the tree level, due to exchange of sneutrinos or 
up- type squarks, as depicted in Figure 1. We need also to include the SM contribution to 



a) 



b) 



FIG. 1: Diagrams for B s — > fj,~ transition within R-parity violating supersymmetric models to 
the lowest order in perturbation theory, (a) due to exchange of sneutrinos, (b) due to exchange of 
up-type squarks. The direction of the sneutrino propagator depends on the helicities of the quark 
and lepton states, in other words whether we have Pl or we have Pr operator at the interaction 
vertex. 



B s — > recall that we are interested in destructive interference of the SM and NP 

amplitudes. Thus, the relevant low-energy \AB\ = 1 effective Hamiltonian would have the 
following form: 



H, 



AB=1 
eff 



TjSM i TTV TTU 

He ft + Hpff + Hp, 



l eff 



eff 



l eff 



Here 
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111 



ttSM 
H eff 



—AGi 



a 



where 



12| 



y/2 2'Ksm 2 9 w 
Y o{xt) = — 



2 n ( V tt V ts) VYYo(x t ) bfP L 8 Vn v P L \l + h.C. 



(10) 



(11) 



x t 



+ 



3x t 



1 - x t ' (1 - x t ) 2 



In x t 



x t = mt/Myy, and r/y is the factor that accounts for the QCD corrections to Yo(x t ). 

Two other terms in Eq. (TlOT) are derived by integrating out the sneutrino and squark heavy 

degrees of freedom. This yields 



H. 



eff 



IPlS p 

mi. 



A ^ Aj22 bP R S fiP L fi + h.C. 



m- 



(12) 



TTU 

H eff 



h^RS WuPLll + h.c. 



(13) 



where z/j , u fc are respectively the lightest sneutrino and the lightest "left" up-type squark 
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states 1 . For (nearly) degenerate sneutrino and/or squark masses, one should replace in 
Eq. ( I12p and/or Eq. f|T3|) the lightest sparticle masses by universal sneutrino and/or squark 
masses, wip^ — > m i > L , rriu kL ~~ >* m u L , as we U as sum over indices i and/or k. 
Also, for the sake of clarity of our analysis, we will follow ref. |j] and assume 

K23 = Km (14) 

in our further calculations. 

Using Eq.'s (HU]) - ( TT3|) as well the simplifying assumption (THj) . one may present the 
B s — > [i + fi~ transition amplitude in the following form: 

A(B S -> M V) = A SM (B S -)■ fj+fi-) + A°{B 8 -> M + /0 + ^"(5, -)■ M V) (15) 



where 



A^(5 S /iV) = -(^Vi^f//I^) = afB Z" (V t tV ts )rj Y Y (x t ) S(p_) 7 6f (P+) (16) 

V 2 sin (7 VK 



A\B S /iV) = -<MVI^//|S.> = ^f^M M(p _) 75 , (p+ ) (18) 

where u(p_) and u(p+) are the bi-spinor wave functions of the leptonic states. (Subsequently, 
p + and p_ are the momenta of fi + and In deriving ffl6|) - ffTHj) we used the following 
parametrization of the hadronic matrix elements: 

(0\IYPls\B s ) = ^p B (0\bj»P R s\B s ) = ^^Pb 



<P\bP R 8\B.) = % (0\bP L s\B s ) ~ lfB ° 



2 \y ////, J ' 2 \_ ////, 

where /# s is the 5 S meson decay constant, and pb is the B s 4- momentum. 

We want to stress that all the three parts of the amplitude have the same structure. They 
all contain the same pseudoscalar bi-spinor bilinear form multiplied by some factor. In what 
follows, both A U (B S — > /i + /U~) and A U (B S — > / u + / u _ ) may interfere with the SM amplitude. 



1 It is assumed that squark mass eigenstates do not differ significantly from the "left" and "right" squark 
states. This is known to be the case for the most of SUSY scenarios with the squark masses much greater 
than 100 GeV. 
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In other words, B s — > fi + fi~ amplitude may have a sign opposite to that of the SM both due 
to the contribution of the sneutrino-mediated diagrams, and due to the contribution of the 
squark-mediated diagram 2 . 

Calculation of the decay branching ratio using ( TT5"j) - (ITSj) is straightforward and yields 



B(B S -> fi+fi- 



tbM b JI 
8vr 



4m 2 
Ml 



G F a m M ,., , 

V2 7rsm z 9 W 



+ 



2m~. m b 



+ 



^2fc2^2fc3 m H 



Am 



(19) 



where tb s is the average lifetime of the neutral B s meson. 



We use during the numerical analysis tb s = 1.509 ps, [131], fs s = 0.225 GeV [l4|, V t lV ts 



0.405 



15| (as mentioned above, we neglect the small phase of this CKM product), Mw 



a 



80.4 GeV, sm 2 9 w = 0.231, G F = 1.166 x 10" 5 GeV~ 
top quark mass we use m v t ole = 173.2 GeV [171 ]. which yields for t 



a(M z ) = 1/128 



16|. For the 



mass mt(m t ) = 163.2 GeV [l|; ?]y = 1.012 for xt 



/■?): 



re MS, QCD renormalized 



We neglect the uncertainties in the values of the input parameters specified above. Those 
are known to alter the predictions for B(B S — > by about 10% [l|, ]2|. This uncertainty 

in B(B S — > fi + fi~) is much less than the one in the experimental value of the branching ratio 
and the one in our results due to destructive interference of different NP amplitudes (see 
the discussion at the end of the paper) . 

We choose mp i > 100 GeV and mj^ > 500 GeV. The squark masses below 500 GeV 
are strongly disfavored by the LHC data (see [16|, [2l| and references therein). To our best 
mowledge, however, no such strong constraints on sneutrino masses has been derived so far 



3. 



Also, following the common approach, we will assume only one non-vanishing RPV cou- 
pling product at a time, or alternatively only one of the NP amplitudes in ( !T4l) to be 
non- vanishing at a time. 

We consider first an idealized scenario with zero uncertainties in the experimental and 
theoretical values of B s — > fi + fi~ branching ratio and perfect coincidence of the Standard 



2 If we give up the simplifying assumption HI]), A V (B S — > will also contain a term with a scalar bi- 

spinor bilinear form. This term however won't interfere with the other terms of the transition amplitude, 
so it does not play any essential role in our analysis. 
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Model prediction with the experimental data. In such a scenario (if assuming real RPV 
couplings), non- vanishing New Physics contribution to B s — > may occur if only the 

transition amplitude has a sign opposite to that of the Standard Model. Following the 
approach of one non-vanishing coupling product at a time, we choose first A 2A , 2 A 2fc3 = or 
equivalently A U (B S — > /U + /i~) = 0. Then the transition amplitude flips the sign if 

A*(B S = -2A SM (B S -> ^V) 

Using Eq.'s (TTB"j) and (117)1 and the values of the input parameters specified above, one finds 
that this occurs when 

-A ?22 A; 23 = 2.12x10-«( i ^) 2 (20) 

This value of A* 22 A- 23 is several times greater in magnitude that the bound quoted in [9j (as 
no amplitude sign flip or large phase has been considered in Q). Nevertheless, Eq. (|20|) im- 
plies rigorous constraints on this coupling product or alternatively on the sneutrino masses. 
Indeed, Eq. (|2"Uj) implies (— A* 22 A- 23 ) ~ for the lightest sneutrino mass ~ 100 GeV. 

Alternatively, if one desires for this coupling product to be of the same order as the SM 
weak coupling squared (g 2 ~ 0.5), the lightest sneutrino should have a mass ~ 50 TeV. This 



is a manifestation of the so-called flavor problem [22|, |23|: to assure that tree level flavor 
changing neutral currents beyond the SM do not conflict with the experimental data, either 
the relevant (effective) couplings should be unnaturally small or the New Physics mass scale 
should be enormously large. Solving the flavor problem goes beyond the scope of the present 
paper. Instead we will simply assume further that A* 22 A- 23 = 0, or A U (B S — > fi + fi~) vanishes. 
Instead we will be concentrating on the contribution of the squark-mediated diagram only 
(Fig. [1] (b)). As mentioned above, possible effects of inference of different NP amplitudes 
will be discussed at the end of the paper. 

If assuming A U (B S — > fi + fi~) = 0, the transition amplitude flips the sign when 

A*(B S -> //+/0 = -2A SM (B S -+ /iV) 

Using Eq.'s ffTBT) and (TTHT) and the values of the input parameters specified above, one finds 
that this occurs when 

- A'oloAU = 6.88 x 10~ 3 u -^— (21) 
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Eq. (l2Tj) implies rather weak constraints on the couplings A' 2fc2 and A 2fc3 . If assuming no 
hierarchy in the values of X' 2k2 and A' 2fc3 , one gets A' 2fc2 ~ 0.085 and X 2k3 ~ 0.085 for ~ 
500 GeV. Thus, moderately small values of A' 2fc2 and A' 2fc3 are still allowed for mu kL ~ 
500 GeV. Furthermore, choosing the lightest "left" up-type squark mass to be heavier (say 
1 TeV of few TeV) would yield larger values for \ 2k2 and \' 2k3 (and for their product) to be 
allowed. 

At first glance this result is not surprising, as the contribution of the diagram with a 
squark exchange in Fig. [TJ (b) is helicity suppressed (like the SM contribution), as can be 
seen e.g. from Eq. ( Tl8l) . We want to stress however that this is a rather non-trivial result, in 
a sense that one should consider the possibility of B s — > n + fi~ amplitude sign flip to derive 
it. If instead one assumes that the sign of the transition amplitude is the same as within 
the SM, so that the NP contribution is solely due to the uncertainty in the experimental 
value of B(B S — > fi + ft~), the constraints on the coupling product X 2k2 \' 2k3 are significantly 
stronger. To illustrate this, we will consider a realistic scenario now: we will demand that 
our predictions for B s — > /i + branching ratio fall in the experimentally allowed interval. 

In order to do this, one should take into account that the experimentally measured 
branching ratio of B s — > fi + ^~ decay is the time integrated branching ratio (it is us ually 
denoted B(B S — > yU + /i~)), which is related to the "theoretical" branching ratio as 0, [3, 25]: 



B(B S -> fM+iM-) = L\j} \ B(B S -> fi+fi-) (22) 

Here 

y s = — ^ = 0.088 ± 0.014 (23) 

zl s 

where Ar s is the width difference in B s — B s mixing, and T s is the average width of the 
neutral B s meson. The expression for A^ T in terms of Wilson coefficients of the low-energy 
effective operators may be found in 2|. For the considered case of real NP couplings and 
under simplifying assumption ([14"]) . one can show after doing some algebra that A^ T = 1. 
Then, using ([2]) and ( 12"3"|) one would get 

(1 - Vs )B{B s fi+fi-) = (2.9±\i) x 10- 9 (24) 

Demanding for the theoretical predictions to fall into the interval specified in the r.h.s. 
of Eq. (|24j) . one gets the following allowed intervals for the coupling product A 2A , 2 A 2fc3 : 

/ 77l~ \ / Tfl~ \ ^ 

- 49X1 °" 500^7 S-^3<9.6X!0- ^ (25) 



and 



5 ' 92 x 10 " [A) £ " A - A - * 7 ' 37 x 10 " [A) (26) 

The first interval (given by ([25]) ) is derived, when the _B S — >■ transition amplitude has 

the same sign as that of the Standard Model. The New Physics contribution is due to the 
uncertainty in the experimental value of B(B S — > /i + /i~). The allowed interval for A 2 £ 2 A 2fc3 
in this case is in a reasonable agreement with that quoted in ref. {9]. The second interval 
(given by ( 1261) ) is derived when the transition amplitude has a sign opposite to that of the 
Standard Model. In that case the allowed values of — A 2fc2 A 2fe3 are g rea t er by an order of 
magnitude. As discussed above, this implies weaker constraints on the allowed region of the 
NP parameter space. 

Notice also that for B s — > ji + n~ amplitude to flip the sign, the coupling product A 2A . 2 A 2fc3 
must be negative (as it follows from Eq. ( 1261) ). Contrary to this, within the other interval 
(given by ( 125]) ). the sign of A 2 * fc2 A 2A . 3 is arbitrary. 

In principle, one may conduct the similar analysis for the contribution of the sneutrino- 
mediated diagrams on Fig. [T] (a). Assuming now that the squark-mediated diagram in 
Fig. [D(b) has a vanishing contribution to B s — > /U + /i~, one will get in this case two different 
intervals for the coupling product A* 22 A^ 23 (that originate in the same way as ( |25]) and ( [26]) 
for A 2fc2 A 2fc3 ). We leave this for a reader as an exercise to do. 

In a realistic scenario neither of the diagrams in Fig. [1] may have a non-vanishing con- 
tribution to B s — > fi + fi~. In addition, one should also take into account the impact of the 
R-conserving sector of the theory on the B s — > fi + n~ transition amplitude as well |6|]. Thus, 
in a realistic scenario one has different sources of a NP contribution to B s — > that 
may in general interfere both constructively and destructively. 

If different NP amplitudes interfere constructively, the coupling product A 2fc2 A 2fc3 may 
also acquire the values between the two intervals given by ( 125]) and ( 126]) . (That is to say, 
the B s — > amplitude sign flip may be only in part due to the contribution of the 

squark-mediated diagram in Fig. [1] (b), it may also be in part due to other New Physics 
effects.) In other words, one should replace ( 125]) and ( 126]) by 

- 4.9 x 1(T 4 ( ^ ) < -A' 2 ; 2 A' 3 < 7.37 x 1(T 3 ( ^ ) (27) 
V500 Ge^y ~ 2fc2 2fc3- \500 GeVJ V ; 

Of course, the NP amplitudes may also interfere destructively. In that case the bounds 
on ^2fc2^2fc3 gi ven by Eq. (127]) may somehow be distorted (they may become weaker). Yet, if 
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FIG. 2: Upper bound on I ^2fc2^2fc3 1 as a function of the NP amplitude phase <3?jvp for rriu^ = 
500 GeV (solid red line), rriu kL = 750 GeV (dashed-dotted green line), rriu^ = 1 TeV (dashed blue 
line). 

there is no fine-tuning or exact cancelation of the contributions of different NP amplitudes, 
it is very unlikely that this distortion alter the bounds on A 2fc2 A 2fc3 , say, by an order of mag- 
nitude. Thus, one may always use ( 12 7 jl to get an insight how large (in order of magnitude) 
the coupling product A 2fc2 A 2fe3 is still allowed to be. 

So far we were assuming that the R-parity violating couplings are real (or have small 
enough phases so that they may be discarded). Yet, our analysis may be extended also to 
the case when these couplings have large phases. Demanding again that our predictions fall 
into the experimentally allowed interval (and assuming again A U (B S — > /i + /U~) to be the only 
non- vanishing NP amplitude), one may derive an upper bound on the absolute value of the 
coupling product A' 2 * fc2 A 2fc3 as a function of the NP amplitude phase $np = arg (A 2 * fc2 A 2fc3 ), 
as depicted in Fig. |2j As one can see from Fig. [2} the bound on |A' 2 * fc2 A 2A , 3 | becomes weaker 
as &np get larger, and it is the weakest when $>np — >■ tt- As mentioned above, this result 
could also be inferred from the analysis of Eq's ([3]) - ([7]). Thus, the most general bound 
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on I^2fc2^2fc3l (both in the case when this product is real and in the case this product has a 
phase) would be 



In conclusion, we have re- visited the problem of New Physics contribution to B s — > /i + yU~ 
decay in light of the recent experimental measurement of this decay branching ratio by the 
LHCb collaboration. We have examined R-parity violating supersymmetric models as a 
primary example, and argued that there is still a room for a significant NP contribution, as 
the transition amplitude still may have a sign opposite to that of the Standard Model or 
alternatively may get a large phase. We have found that if taking into account the effect 
of B s — » amplitude possible sign flip as compared to that of the SM (or possible 

large phase), the bounds imposed on the RPV coupling products that account for the 
B s — > fi + 'yT transition may be weaker by an order of magnitude than if the effect of the 
amplitude sign flip (or possible large phase) is disregarded. We emphasize that a similar 
effect may be observed also within other SM extensions. So considering within other New 
Physics models the possibility for the B s — > y + y~ transition amplitude to have a sign 
opposite to that of the SM or have a large phase is strongly encouraged. 

The author is grateful to Alexey A. Petrov for stimulating discussions and valuable com- 
ments. 
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